Transport properties of n-type ultrananocrystalline diamond films 



O 

o 

O 

Q 

(N 



I. S. Beloborodov 1 ' 2 , P. Zapol 1,3,4 , D. M. Gruen 1 , and L. A. Curtiss 1,3,4 
1 Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 
2 Department of Physics, University of Chicago, Chicago, Illinois 60637, USA 
3 Center for Nanoscale Materials, Argonne National Laboratory, Argonne, Illinois 60439, USA 
4 Chemistry Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 

(Dated: February 6, 2008) 

We investigate transport properties of ultrananocrystalline diamond films for a broad range of 
temperatures. Addition of nitrogen during plasma-assisted growth increases the conductivity of 
ultrananocrystalline diamond films by several orders of magnitude. We show that films produced 
at low concentration of nitrogen in the plasma are very resistive and electron transport occurs via 
a variable range hopping mechanism while in films produced at high nitrogen concentration the 
electron states become delocalized and the transport properties of ultrananocrystalline diamond 
films can be described using the Boltzmann formalism. We discuss the critical concentration of 
carriers at which the metal to insulator transition in ultrananocrystalline diamond films occurs and 
compare our results with available experimental data. 

PACS numbers: 81.05.Uw, 73.22-f, 73.63.Bd, 72.10.-d 
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I. INTRODUCTION 

Ultrananocrystalline diamond (UNCD) films have 
many potential applications in high-temperature elec- 
tronics because of the highest n-type conductivity and 
carrier concentration demonstrated for phase-pure dia- 
mond thin films. 1,2,3 Experimentally, n-type conductivity 
of UNCD can be tuned within several orders of magni- 
tude by nitrogen addition to the feed gas during plasma- 
assisted growth. Without the nitrogen addition to the 
plasma, the resulting films are quite resistive, but addi- 
tion of more than about 8% of nitrogen leads to a metal- 
lic behavior. Interestingly, the grain size increases from 
4 nm to 16 nm, as well as the grain boundary width 
from about 0.5 nm to 2.2 nm with N concentration in- 
crease from 0% to 20%, but the overall fraction of dia- 
mond phase remains consistently very high.— At the same 
time, nitrogen concentration in the films remains nearly 
the same, suggesting that change in the carrier concen- 
tration due to higher number of nitrogen donors cannot 
be invoked to explain the increase of conductivity. It was 
shown by careful Hall measurements- that the conductiv- 
ity is n-type, and that carrier concentrations in the films 
increase much more than carrier mobilities with increase 
in feed gas nitrogen concentration. Therefore, these stud- 
ies suggested that there is a metal-insulator transition in 
UNCD upon nitrogen addition, and that this transition is 
triggered by changes in UNCD grain and grain boundary 
sizes. 

Theoretical tight-binding studies^ of UNCD with- 
out nitrogen suggested a hopping transport mechanism, 
where the donor sites are located in the grain bound- 
aries and the bulk of the grains is electrically inactive. 
Here high- angle high-energy grain boundaries were con- 
sidered to be a model of general grain boundaries of 
UNCD. The donor states were calculated to be carbon tt- 
states and dangling bond states in the grain boundaries. 
The hopping conductivity of UNCD was estimated to be 
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FIG. 1: Calculated DOS for diamond with E13 grain bound- 
ary and without (left) or with (right) nitrogen (relative to the 
Fermi level of undoped system). Occupied states are shaded. 
Sketch of electron levels in the presence of nitrogen is shown 
on the right. The data are taken from Ref. 0- 



10 -6 (f2cm) -1 , 3 to 7 orders of magnitude higher than in 
CVD diamond. 

Self-consistent charge density functional based tight- 
binding studies 7 of nitrogen sites in the bulk and in the 
grain boundaries found energies, geometries and elec- 
tronic structure of nitrogen in UNCD at low and high 
nitrogen concentrations. First, it was found that nitro- 
gen is energetically preferred in the grain boundary rel- 
ative to the bulk diamond sites in the grains. Second, 
nitrogen was found not to be a shallow donor, being pref- 
erentially in non-doping configurations. Third, electrons 
from high-lying donor states formed upon nitrogen addi- 
tion transfer down in energy to the Fermi level, located 
in the 7r-state region. This causes an upward shift of 
the Fermi level with increase in nitrogen concentration. 
Calculated DOS for diamond grain boundaries with and 
without nitrogen are shown in fig. [T] along with a scheme 
of electron levels in the latter case. In addition, the con- 
nectivity of the 7r-bonds in the grain boundaries improves 
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FIG. 2: Sketch of ultrananocrystalline diamond film with 
a schematic electron path. The grey regions are insulating 
grains with a and d being the characteristic grain size and 
distance between the grains respectively. 



upon nitrogen addition. This Fermi level shift and im- 
proved connectivity were proposed to be responsible for 
increase in hopping conductivity. 

However, there is still no satisfactory mesoscopic de- 
scription of UNCD transport properties. No explanation 
was provided previously for the transition from hopping 
to semimetallic regime with increase in N concentration. 
In this paper we consider theoretically transport proper- 
ties of UNCD films in the temperature range 



10K < T < 700K, 



(1) 



and investigate UNCD conductivity within approxima- 
tions developed for an array of metallic and semiconduct- 
ing grains^ in both Ohmic and non-Ohmic regimes. 



II. INSULATING REGIME 

For low concentrations of nitrogen (< 1% in the feed 
gas), the resistance of UNCD films per square is R ~ 
10 - 100 MVL > R c , where R c = is the quantum 

resistance. In this case, UNCD films are insulators. 

The characteristic grain sizes, see Fig. [5J of UNCD 
films are in the range 

Anm < a < 16nm (2) 

and characteristic width of the grain boundaries is 

0.5nm < d < l.lnm. (3) 

The charging energies E c — e 2 / 'nd, with k being the di- 
electric constant, corresponding to these grain boundary 
widths are 



- 10 3 K < E c < - 4 x 10 4 K. 

K K 



(4) 



Assuming that the dielectric constant for UNCD films is 
k ~ 3 — 6 from Eq. follows that the temperature T 
satisfies the inequality 

T < E c . (5) 
In this temperature range the activation conductivity, 

a ~ exp(-£: c /T) (6) 



(tunneling between the nearest neighbors), is suppressed 
and electron transport occurs through tunneling via sev- 
eral grains. 

Below we assume that the concentration of donors in 
the grain boundaries, n^, satisfies the condition 

«^ 3 «1, (7) 

where £ is the localization length. Using the fact that 

f ~ d ~ 1 nm, (8) 

we obtain the following condition for concentration of 
donors 



n d < 10 cm" 



(9) 



We first consider phonon assisted (Ohmic) transport in 
UNCD films and later discuss the non-Ohmic regime. 



A. Ohmic transport 



Following the Mott-Efros-Shklovskii theor; 



.10.11 



we 



write the hopping conductivity in the following form 



=2 \ 



(10) 



where £ is the localization length (£ is of the order of a 
boundary width, d). Equation () 1 0(1 is valid for 3D sam- 
ples and thick films with the sample thickness L larger 
than the optimal hopping length r opt defined below. This 
is a typical case for UNCD films, which are grown to 
thicknesses of fim^. Minimizing Eq. (JTUJ) with respect to 
the hopping distance, r, we obtain the optimal hopping 
length 



' opt 




(11) 



where To = e 2 /k£ is the characteristic energy scale. From 
Eq. (fTTj) follows that the hopping length r opt (T) increases 
with decreasing temperature. Using the fact that the 
hopping length 

r opt (T) = N(T)a, (12) 
where N(T) is the number of grains in one hop, we obtain 



(13) 



For temperatures T < 300K and T - 40000if one ob- 
tains 



N 



(14) 



With increase in temperature the number of grains N in 
one hop will decrease, leading to the hopping distance 
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r pt(T) of the order of the grain size. Also, increase in 
nitrogen concentration is shown to lead to increase in 
the grain size. With increase in the grain size a (and 
proportional increase in the grain boundary width d) N 
will decrease as ^fa. Substituting Eq. (fTTj) into Eq. (flQ|) 
we finally obtain Efros-Shklovskii typei^ of conductivity 
behavior 



exp 



(-y/n/f). 



(15) 



Equation (I15|) is valid for temperatures T < To- At tem- 
peratures T ~ To the hopping distance r op t(T) becomes 
of the order of a grain size and Eq. fTO]) does not work 
anymore. In this temperature range the conductivity 
obeys the Arrhenius law (J6j) . 

We would like to discuss the possibility of the observa- 
tion of the Mott law in UNCD films. In semiconductors, 
the Efros-Shklovskii law may turn to the Mott behavior 
with the increase of temperature. This happens when 
the typical electron energy e involved in a hopping pro- 
cess becomes larger than the width of the Coulomb gap 
A c , i.e. when it falls into the flat region of the density 
of states where Mott behavior is expected. To estimate 
the width of the Coulomb gap A c - one compares the ES 
expression for the density of states 



iD-l 



(16) 



where D is the dimensionality of a sample, with the bare 
density of states vq, i.e. the DOS in the absences of the 
long-range part of the Coulomb interactions. Using the 
condition v(A c ) ~ vq we obtain 



A, = 



u e 



2D 



Inserting the value for the bare DOS, v 
Eq. (fT7|) we finally obtain 



(17) 



l/E c £ D , into 



(18) 



Equation (jTSJ) means that there is no flat region in the 
density of ground states and, thus, the Mott regime is 
difficult to observe in UNCD films. At first glance, this 
conclusion contradicts numerically obtained results us- 
ing tight-binding methods^ since non-zero densities of 
states at the Fermi level was obtained as a result of these 
calculations. However, there are relatively few states in 
the numerical studies because of the number of atom lim- 
itations on the size of the unit cells, which leads to uncer- 
tainties in the peak positions larger than typical values 
of Coulomb gap of several 10K. In addition, Coulomb 
gap existence is based on excitonic effects, which are not 
explicitly taken into account in the tight-binding simula- 
tions. In other words, tight-binding simulations are not 
designed to accurately reproduce the functional form of 
DOS in the immediate vicinity (within several tens of 
Kelvins) of the Fermi level. 



B. Non-Ohmic transport 

Here we consider non-Ohmic conductivity in the tem- 
perature range in which the Ohmic conductivity is given 
by Eq. (fl3|) . In the presence of electric field, E, the hop- 
ping conductivity is given by the following expression 



a ~ exp 



ktT 



eEr 
~T~ 



(19) 



1. Weak electric field 

In a weak electric field, eE^ < T, the last term in 
the r.h.s. of Eq. (fT9| is not important and the optimal 
hopping length is given by Eq. (fTT|) . Using Eqs. (fT9|) and 
(|11[) one obtains for the hopping conductivity in a weak 
electric field the following expression 



exp exp [(eE£/T) ^%Jt 



(20) 



From Eq. (|20|) follows that the Coulomb term becomes 
of the order of the electric field term at fields 



E* 



T/e£. 



(21) 



2. Strong electric field 

For electric fields E > E* , where the field E* is defined 
by Eq. (|2ip . the first term in the r.h.s. of Eq. (fT9|) is not 
important. The optimal hopping length in this regime is 
temperature independent and is given by the following 
expression 



' opt 



(E) = ZVEtjE, 



(22) 



where E = e/< 2 = T /e£ (we note that E > E*). 
Using Eq. (f22|) one obtains for the hopping conductivity 
in a strong electric fielcU^ 



exp 



■VEo/E 



(23) 



We note that: i) for very strong electric fields, E ~ Eq, 
where the field Eq is defined below Eq. (|22|) . the hop- 
ping picture does not work anymore since the hopping 
length becomes of the order of the localization length, 
fopt{Eo) ~ £; ii) The current density j can be expressed 
in term of conductivity as follows j — aE, where the 
conductivity a is given by Eqs. (|20p or (f2"3")) depending 
on the strength of the electric field E. 



III. METALLIC REGIME 

The transition from the hopping regime to degener- 
ate semiconductor regime in UNCD is associated with 
a strong increase in the number of carriers, whereas the 
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mobility of carriers is not strongly affected. This was 
demonstrated in Hall measurements^. Also, higher ni- 
trogen concentration is related to increase in grain size 
and grain boundary width. For relatively large concen- 
trations of nitrogen in the plasma, (3> 5%), the UNCD 
film resistance becomes less than the quantum resistance, 
R ~ 100 fl <C R c . As a result, the electron states become 
delocalized and the transport properties of UNCD films 
can be described using the Boltzmann equation. 

There are several scattering mechanisms which con- 
tribute to the resistance, R, of a sample 14 : 1) scattering 
of electrons on static impurities, R e i, 2) electron-phonon 
scattering, i? e _ p ; l , and 3) electron-electron scattering, 
R e - e . At high temperatures the interference between 
different mechanisms can be neglected and the total re- 
sistance consists of a sum of resistances 



R — R P 



where i? P , = ar„- , R P 



R P 



-ph 



Re 



(24) 



'*T-life, Re 



le-ph — <-*i e —ph> -"-e-e 

Here a — m*/(e 2 n e ) with m* being the effective mass of 
an electron and r e j, r e _ p /j and r e _ e are the characteristic 
scattering times due to scattering of electrons on static 
impurities, phonons and electrons respectively. We would 
like to comment on the validity of Eq. (|24p. Writing 
Eq. I|24p we assume that different scattering mechanisms 
are independent. This is indeed the case at not very low 
temperatures when the scattering lengths L e ^ p h = u/T, 
with u being the phonon speed and L e - e = vf/T are 
smaller than the clastic mean fre6 path., / — VF^~ei' 

Us- 
ing the condition max(L e _ p ^, L e -e) ~ I one can find the 
characteristic temperature T* = . For temperatures 
T > T* the interference among different mechanisms is 
not important and Eq. (I24| is valid. Below we discuss 
different contributions to the resistance of UNCD films. 

Assuming that the mean free path of an electron is 
of the order of the grain boundary width d for elastic 
scattering time one obtains r" 1 ~ jj^/d, where Vp is 
the Fermi velocity. The electron-phonon scattering time 
T e -ph is given by the following expression™ 



e—ph 




T > LO D 
T < LJ D , 



(25) 



where lod is the Debye temperature (for UNCD films 
lod ~ 1000 if). For the temperature range considered 
here, only the T < cup regime is relevant. For electron- 
electron scattering time r e _ e one finds rjl e = T 2 //i, 
where \i is the chemical potential. Using Eqs. (|24|) and 
(|2"B")) for temperatures T < top, we obtain the following 
expression for the total resistance R of a sample 

R = {m*/e 2 n e ) [v F /d + T 5 /lo 4 d + T 2 /fi}. (26) 

For temperatures T < top the last two terms in the r.h.s. 
of Eq. I|26p are small therefore the resistance R is almost 
temperature independent. 



IV. CROSSOVER FROM INSULATING TO 
METALLIC REGIME 

In this section we discuss crossover from insulating 
(hopping) to semimetallic regime in UNCD films. To 
find a critical concentration of carriers n c at which the 
crossover from insulating to semimetallic regime occurs 
following the Mott's approach 1 ^, we compare two char- 
acteristic lengths: the screening length r t 



° 2 (£3„)l/6 

and the localization length £. In the metallic regime the 
screening length r e is small and the localization length £ 
is large while in the insulating regime the situation is the 
opposite. The crossover region between the two regimes 
can be defined as r e ~ £. Using the above condition we 
obtain the following expression for the critical concentra- 
tion of carriers 



1 



(40 s 



(27) 



Since the localization length £ is of the order of the 
distance between the grains, £ ~ d, see Eq. ((5J, from 
Eq. (|27|) follows that the larger the distance between the 
grains the smaller the critical concentration n c is needed 
to achieve the metal to insulator crossover. For a realistic 
d ~ 1 nm we obtain the following estimate for the critical 
concentration n c ~ (10 19 — 10 20 ) cm~ 3 , which appears to 
be plausible. 2 



V. COMPARISON WITH EXPERIMENT 

In this section we compare our results for resistivity 
of UNCD films with available experimental data from 
Refs. and y. These data are shown in Fig. [3l One can 
distinguish two qualitatively different regimes for con- 
ductivity behavior: i) metallic, a > 117 _1 cm _1 , and ii) 
semiconducting, a < li7 _1 cm _1 . 

To compare our results for resistivity R in the metal- 
lic regime, Eq. (|2"6"|) . with experimental values we have 
taken reasonable estimates for the critical carrier con- 
centration of n e ~ n c ~ 10 20 cm~ 3 and vp of 10 8 cm/s. 
From Eq. (f26|) , we obtain the specific resistivity of UNCD 
to be R ~ 0.1 fkm. We note that Eq. (|2"6"|) was estimated 
to be valid above temperatures T* = 10 — 50 K. 

In the semiconducting regime the conductivity behav- 
ior of UNCD films is described by Eq. (fT5]) . Since the en- 
ergy scale To in Eq. (fT5|) is inversely proportional to the 
localization length £ one can see that at fixed tempera- 
ture T the conductivity a is very sensitive to the value 
of localization length. This fact is in perfect agreement 
with experimental data in Fig. [3J 

The transport properties of UNCD films have been dis- 
cussed in several publications^^. Our treatment does 
not contradict the proposed DOS behavior 3 . In fact, 
electron densities in this work are consistent with our 
treatment. The main difference between previous con- 
siderations and our work is based on our description of 
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FIG. 3: Dependence of the conductivity of UNCD films on 
temperature. Different curves correspond to different nitro- 
gen concentrations in the feed gas during film growth. Repro- 
duced from Ref. @ 



experimental data using two regimes i) extended elec- 
tronic states, Sec. IIIII and ii) localized electronic states, 



Sec. |TTJ All previous interpretations of experimental data 
were based on electron hopping theory only. 



VI. CONCLUSIONS 

We studied transport properties of UNCD films as a 
function of temperature and concentration of nitrogen 
during plasma-assisted growth. We show that for UNCD 
grown in low concentration of nitrogen (< 1% in the feed 
gas) the electron transport occurs via a variable range 
hopping mechanism while for films grown at high nitro- 
gen concentration the electron states become delocalized 
and the transport properties of ultrananocrystalline dia- 
mond films can be described using the Boltzmann equa- 
tion. Our results are in agreement with available experi- 
mental data. 
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